A complete map of specificity encoding for a partially fuzzy protein interaction.

Thousands of human proteins function by binding short linear motifs embedded in intrinsically disordered regions. How affinity and specificity are encoded in these binding domains and the motifs themselves is not well understood. The evolvability of binding specificity - how rapidly and extensively it can change upon mutation - is also largely unexplored, as is the contribution of 'fuzzy' dynamic residues to affinity and specificity in protein-protein interactions. Here we report the first complete map of specificity encoding for a globular protein domain. Quantifying >200,000 energetic interactions between a PDZ domain and its ligand identifies 20 major energetically coupled pairs of sites that control specificity. These are organized into six modules, with most mutations in each module reprogramming specificity for a single position in the ligand. Nine of the major energetic couplings controlling specificity are between structural contacts and 11 have an allosteric mechanism of action. The dynamic tail of the ligand is more robust to mutation than the structured residues but contributes additively to binding affinity and communicates with structured residues to enable changes in specificity. Our results quantify the binding specificities of >1,800 globular proteins to reveal how specificity is encoded and provide a direct comparison of the encoding of affinity and specificity in structured and dynamic molecular recognition.

Evolution of short linear motifs and disordered proteins Topic: yeast as model system to study evolution.

Evolutionary preservation of protein structure had a major influence on the field of molecular evolution: changes in individual amino acids that did not disrupt protein folding would either have no effect or subtly change the 'lock' so that it could fit a new 'key'. Homology of individual amino acids could be confidently assigned through sequence alignments, and models of evolution could be tested. This view of molecular evolution excluded large regions of proteins that could not be confidently aligned, such as intrinsically disordered regions (IDRs) that do not fold into stable structures. In the last decade, major progress has been made in understanding the evolution of IDRs, much of it facilitated by new experimental and computational approaches in yeast. Here, we review this progress as well as several still outstanding questions.

Discovering molecular features of intrinsically disordered regions by using evolution for contrastive learning.

A major challenge to the characterization of intrinsically disordered regions (IDRs), which are widespread in the proteome, but relatively poorly understood, is the identification of molecular features that mediate functions of these regions, such as short motifs, amino acid repeats and physicochemical properties. Here, we introduce a proteome-scale feature discovery approach for IDRs. Our approach, which we call "reverse homology", exploits the principle that important functional features are conserved over evolution. We use this as a contrastive learning signal for deep learning: given a set of homologous IDRs, the neural network has to correctly choose a held-out homolog from another set of IDRs sampled randomly from the proteome. We pair reverse homology with a simple architecture and standard interpretation techniques, and show that the network learns conserved features of IDRs that can be interpreted as motifs, repeats, or bulk features like charge or amino acid propensities. We also show that our model can be used to produce visualizations of what residues and regions are most important to IDR function, generating hypotheses for uncharacterized IDRs. Our results suggest that feature discovery using unsupervised neural networks is a promising avenue to gain systematic insight into poorly understood protein sequences.

Identifying molecular features that are associated with biological function of intrinsically disordered protein regions.

In previous work, we showed that intrinsically disordered regions (IDRs) of proteins contain sequence-distributed molecular features that are conserved over evolution, despite little sequence similarity that can be detected in alignments (Zarin et al., 2019). Here, we aim to use these molecular features to predict specific biological functions for individual IDRs and identify the molecular features within them that are associated with these functions. We find that the predictable functions are diverse. Examining the associated molecular features, we note some that are consistent with previous reports and identify others that were previously unknown. We experimentally confirm that elevated isoelectric point and hydrophobicity, features that are positively associated with mitochondrial localization, are necessary for mitochondrial targeting function. Remarkably, increasing isoelectric point in a synthetic IDR restores weak mitochondrial targeting. We believe feature analysis represents a new systematic approach to understand how biological functions of IDRs are specified by their protein sequences.

Whence Blobs? Phylogenetics of functional protein condensates.

What do we know about the molecular evolution of functional protein condensation? The capacity of proteins to form biomolecular condensates (compact, protein-rich states, not bound by membranes, but still separated from the rest of the contents of the cell) appears in many cases to be bestowed by weak, transient interactions within one or between proteins. Natural selection is expected to remove or fix amino acid changes, insertions or deletions that preserve and change this condensation capacity when doing so is beneficial to the cell. A few recent studies have begun to explore this frontier of phylogenetics at the intersection of biophysics and cell biology.

Condensation of Ded1p Promotes a Translational Switch from Housekeeping to Stress Protein Production.

Cells sense elevated temperatures and mount an adaptive heat shock response that involves changes in gene expression, but the underlying mechanisms, particularly on the level of translation, remain unknown. Here we report that, in budding yeast, the essential translation initiation factor Ded1p undergoes heat-induced phase separation into gel-like condensates. Using ribosome profiling and an in vitro translation assay, we reveal that condensate formation inactivates Ded1p and represses translation of housekeeping mRNAs while promoting translation of stress mRNAs. Testing a variant of Ded1p with altered phase behavior as well as Ded1p homologs from diverse species, we demonstrate that Ded1p condensation is adaptive and fine-tuned to the maximum growth temperature of the respective organism. We conclude that Ded1p condensation is an integral part of an extended heat shock response that selectively represses translation of housekeeping mRNAs to promote survival under conditions of severe heat stress.

Proteome-wide signatures of function in highly diverged intrinsically disordered regions.

Intrinsically disordered regions make up a large part of the proteome, but the sequence-to-function relationship in these regions is poorly understood, in part because the primary amino acid sequences of these regions are poorly conserved in alignments. Here we use an evolutionary approach to detect molecular features that are preserved in the amino acid sequences of orthologous intrinsically disordered regions. We find that most disordered regions contain multiple molecular features that are preserved, and we define these as 'evolutionary signatures' of disordered regions. We demonstrate that intrinsically disordered regions with similar evolutionary signatures can rescue function in vivo, and that groups of intrinsically disordered regions with similar evolutionary signatures are strongly enriched for functional annotations and phenotypes. We propose that evolutionary signatures can be used to predict function for many disordered regions from their amino acid sequences.

YeastSpotter: accurate and parameter-free web segmentation for microscopy images of yeast cells.

SUMMARY: We introduce YeastSpotter, a web application for the segmentation of yeast microscopy images into single cells. YeastSpotter is user-friendly and generalizable, reducing the computational expertise required for this critical preprocessing step in many image analysis pipelines. AVAILABILITY AND IMPLEMENTATION: YeastSpotter is available at http://yeastspotter.csb.utoronto.ca/. Code is available at https://github.com/alexxijielu/yeast_segmentation. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Short linear motifs in intrinsically disordered regions modulate HOG signaling capacity.

BACKGROUND: The effort to characterize intrinsically disordered regions of signaling proteins is rapidly expanding. An important class of disordered interaction modules are ubiquitous and functionally diverse elements known as short linear motifs (SLiMs). RESULTS: To further examine the role of SLiMs in signal transduction, we used a previously devised bioinformatics method to predict evolutionarily conserved SLiMs within a well-characterized pathway in S. cerevisiae. Using a single cell, reporter-based flow cytometry assay in conjunction with a fluorescent reporter driven by a pathway-specific promoter, we quantitatively assessed pathway output via systematic deletions of individual motifs. We found that, when deleted, 34% (10/29) of predicted SLiMs displayed a significant decrease in pathway output, providing evidence that these motifs play a role in signal transduction. Assuming that mutations in SLiMs have quantitative effects on mechanisms of signaling, we show that perturbations of parameters in a previously published stochastic model of HOG signaling could reproduce the quantitative effects of 4 out of 7 mutations in previously unknown SLiMs. CONCLUSIONS: Our study suggests that, even in well-characterized pathways, large numbers of functional elements remain undiscovered, and that challenges remain for application of systems biology models to interpret the effects of mutations in signaling pathways.

Parallel reorganization of protein function in the spindle checkpoint pathway through evolutionary paths in the fitness landscape that appear neutral in laboratory experiments.

Regulatory networks often increase in complexity during evolution through gene duplication and divergence of component proteins. Two models that explain this increase in complexity are: 1) adaptive changes after gene duplication, such as resolution of adaptive conflicts, and 2) non-adaptive processes such as duplication, degeneration and complementation. Both of these models predict complementary changes in the retained duplicates, but they can be distinguished by direct fitness measurements in organisms with short generation times. Previously, it has been observed that repeated duplication of an essential protein in the spindle checkpoint pathway has occurred multiple times over the eukaryotic tree of life, leading to convergent protein domain organization in its duplicates. Here, we replace the paralog pair in S. cerevisiae with a single-copy protein from a species that did not undergo gene duplication. Surprisingly, using quantitative fitness measurements in laboratory conditions stressful for the spindle-checkpoint pathway, we find no evidence that reorganization of protein function after gene duplication is beneficial. We then reconstruct several evolutionary intermediates from the inferred ancestral network to the extant one, and find that, at the resolution of our assay, there exist stepwise mutational paths from the single protein to the divergent pair of extant proteins with no apparent fitness defects. Parallel evolution has been taken as strong evidence for natural selection, but our results suggest that even in these cases, reorganization of protein function after gene duplication may be explained by neutral processes.

Selection maintains signaling function of a highly diverged intrinsically disordered region.

Intrinsically disordered regions (IDRs) are characterized by their lack of stable secondary or tertiary structure and comprise a large part of the eukaryotic proteome. Although these regions play a variety of signaling and regulatory roles, they appear to be rapidly evolving at the primary sequence level. To understand the functional implications of this rapid evolution, we focused on a highly diverged IDR in Saccharomyces cerevisiae that is involved in regulating multiple conserved MAPK pathways. We hypothesized that under stabilizing selection, the functional output of orthologous IDRs could be maintained, such that diverse genotypes could lead to similar function and fitness. Consistent with the stabilizing selection hypothesis, we find that diverged, orthologous IDRs can mostly recapitulate wild-type function and fitness in S. cerevisiae We also find that the electrostatic charge of the IDR is correlated with signaling output and, using phylogenetic comparative methods, find evidence for selection maintaining this quantitative molecular trait despite underlying genotypic divergence.

Turnover of protein phosphorylation evolving under stabilizing selection.

Most proteins are regulated by posttranslational modifications and changes in these modifications contribute to evolutionary changes as well as to human diseases. Phosphorylation of serines, threonines, and tyrosines are the most common modifications identified to date in eukaryotic proteomes. While the mode of action and the function of most phosphorylation sites remain unknown, functional studies have shown that phosphorylation affects protein stability, localization and ability to interact. Two broad modes of action have been described for protein phosphorylation. The first mode corresponds to the canonical and qualitative view whereby single phosphorylation sites act as molecular switches that either turn on or off specific protein functions through direct or allosteric effects. The second mode is more akin to a rheostat than a switch. In this case, a group of phosphorylation sites in a given protein region contributes collectively to the modification of the protein, irrespective of the precise position of individual sites, through an aggregate property. Here we discuss these two types of regulation and examine how they affect the rate and patterns of protein phosphorylation evolution. We describe how the evolution of clusters of phosphorylation sites can be studied under the framework of complex traits evolution and stabilizing selection.

Cyclosporin A enhances neural precursor cell survival in mice through a calcineurin-independent pathway.

Cyclosporin A (CsA) has direct effects on neural stem and progenitor cells (together termed neural precursor cells; NPCs) in the adult central nervous system. Administration of CsA in vitro or in vivo promotes the survival of NPCs and expands the pools of NPCs in mice. Moreover, CsA administration is effective in promoting NPC activation, tissue repair and functional recovery in a mouse model of cortical stroke. The mechanism(s) by which CsA mediates this cell survival effect remains unknown. Herein, we examined both calcineurin-dependent and calcineurin-independent pathways through which CsA might mediate NPC survival. To examine calcineurin-dependent pathways, we utilized FK506 (Tacrolimus), an immunosuppressive molecule that inhibits calcineurin, as well as drugs that inhibit cyclophilin A-mediated activation of calcineurin. To evaluate the calcineurin-independent pathway, we utilized NIM811, a non-immunosuppressive CsA analog that functions independently of calcineurin by blocking mitochondrial permeability transition pore formation. We found that only NIM811 can entirely account for the pro-survival effects of CsA on NPCs. Indeed, blocking signaling pathways downstream of calcineurin activation using nNOS mice did not inhibit CsA-mediated cell survival, which supports the proposal that the effects are calcinuerin-independent. In vivo studies revealed that NIM811 administration mimics the pro-survival effects of CsA on NPCs and promotes functional recovery in a model of cortical stroke, identical to the effects seen with CsA administration. We conclude that CsA mediates its effect on NPC survival through calcineurin-independent inhibition of mitochondrial permeability transition pore formation and suggest that this pathway has potential therapeutic benefits for developing NPC-mediated cell replacement strategies.

Insights into molecular evolution from yeast genomics.

Enabled by comparative genomics, yeasts have increasingly developed into a powerful model system for molecular evolution. Here we survey several areas in which yeast studies have made important contributions, including regulatory evolution, gene duplication and divergence, evolution of gene order and evolution of complexity. In each area we highlight key studies and findings based on techniques ranging from statistical analysis of large datasets to direct laboratory measurements of fitness. Future work will combine traditional evolutionary genetics analysis and experimental evolution with tools from systems biology to yield mechanistic insight into complex phenotypes.